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ABSTRACT 

A general thermodynamic argument shows that multi-component matter in full 
chemical equilibrium, with uniform entropy per baryon, is generally stably stratified. 
This is particularly relevant for neutron stars, in which the effects of entropy are neg- 
ligible compared to those of the equilibrium composition gradient established by weak 
interactions. It can therefore be asserted that, regardless of the uncertainties in the 
equation of state of dense matter, neutron stars are stably stratified. This has im- 
portant, previously discussed consequences for their oscillation modes, magnetic field 
evolution, and internal angular momentum transport. 

Subject headings: convection — dense matter — equation of state — hydrodynamics 
— stars: neutron — stars: oscillations 

1. Introduction 

Whether a fluid is stably stratified or unstable to convective overturn plays a crucial role in its 
dynamics and transport properties. This is particularly true in non-degenerate stars, where heat 
may be transported by convection or by radiation, and which of the two mechanisms is effective 
determines the shape of the temperature profile, and with it much of the structure of the star (see, 
e.g., Kippenhahn & Weigert (1994)). 

In the case of neutron stars, this issue had been largely ignored until it was pointed out that 
the transport of magnetic flux could be strongly affected if the matter was stably stratified (Pethick 
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(1992); Goldreich Sz Reisenegger (1992)). Simple models for the composition and equation of state 
of equilibrium neutron star matter generally show that it is stably stratified due to the presence of a 
chemical gradient, and give more or less consistent estimates for the characteristic (Brunt- Vaisala) 
frequency (Pethick (1992); Reisenegger & Goldreich (1992); Lai (1994); Lee (1995)). 

Stable stratification leads to a set of g-modes (Reisenegger & Goldreich (1992); Lai (1994)) 
which might be excited during binary inspiral (Reisenegger & Goldreich (1994); Lai (1994); Ho 
& Lai (1999)). It also may affect angular momentum transport by inhibiting Ekman circulation 
(Reisenegger (1993, 1995); Abney & Epstein (1996)), with possible observational signatures in 
pulsar glitches (Abney, Epstein, & Olinto (1996)). Its effects on magnetic flux evolution have 
also gained relevance with the interpretation of soft gamma-ray repeaters and anomalous X-ray 
pulsars as ultramagnetized neutron stars powered by magnetic field decay (Thompson & Duncan 
(1995, 1996); Kouveliotou et al. (1998); Gotthelf, Vasisht, & Dotani (1999)). The recent work on r- 
modes destabilized by gravitational wave emission in very rapidly rotating neutron stars (Andersson 
(1998); Lindblom, Owen, & Morsink (1998)) is also paying increasing attention to the effects of 
stable stratification and the relation between the modes of stratified and unstratified stars (e.g., 
Lockitch & Friedman (1999)). 

The previous work has made it plausible that cold neutron star matter is always stably strat- 
ified, but it has only considered specific models for neutron star matter, which are still quite 
uncertain (Lattimer & Prakash (2000)). The present paper shows by a simple thermodynamic ar- 
gument that an isentropic, multi-component fluid in full chemical equilibrium cannot be unstable to 
convection. 2 Since cold neutron stars are zero-entropy objects, this proves the point above, except 
for the remaining, but unlikely possibility of marginal stability (discussed below). 

In the next section, the stage is set by a heuristic derivation of the general condition for stable 
stratification, given previously by other authors, but now phrased in the language and notation to 
be used in the rest of this work. Section 3 contains the thermodynamic proof that is the core of 
the present paper. Implications of this result are discussed in the final section. 



2. Stable stratification 

The stability condition for a static fluid in a gravitational field 3 g = — can be obtained 
by the following heuristic argument (e.g., Cox (1980)). The (Newtonian) condition of hydrostatic 
equilibrium, 

VP = - P V0, (1) 



2 A different argument was recently given by Thompson (2000). 

3 In the case of a rotating star with a simple rotation law, the gravitational field and the corresponding potential 
can be generalized to include centrifugal effects. 
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requires the gravitational potential cf>, the pressure P, and the mass density p in the equilibrium 
configuration of the fluid to be constant on the same set of surfaces. Take a fluid element and 
displace it perpendicularly to these surfaces, on a timescale slow enough to keep it in pressure 
equilibrium with the surrounding matter, but not necessarily slow enough to exchange heat and 
particles with its surroundings or reach chemical equilibrium by internal processes. Therefore, the 
fluid element will generally reach a different density than the surrounding matter, and a buoyancy 
force will either increase or decrease its initial displacement. The buoyancy force per unit mass 
per unit displacement (taken to be positive if it acts so as to restore the fluid element to its initial 
position) is the squared Brunt- Vaisala frequency, 



N = 



dP) eq \dP) P ert 



g\\ (2) 



where the subscripts "eg" and "pert" refer, respectively, to variations in the background equilibrium 
fluid and in the perturbed fluid element. If A^ 2 > 0, the fluid can support oscillations with buoyant 
restoring forces (gravity waves) and characteristic frequency N. If N 2 < 0, the fluid is unstable to 
convective overturn. 

Take the entropy per baryon of the equilibrium fluid to be S(P), and its chemical composition 
to be parameterized by a set of variables Yi(P) (abundance, per baryon, of particle species i = 
1,2, ...,n). The pressure is used to mark the depth in the fluid at which the variables are to be 
evaluated. Assume that the displaced fluid element does not change its internal variables S and 
{Yi}f =1 . The condition for stable stratification, A^ 2 > 0, can then be rewritten in the alternative 
forms 

dp \ f dp 



>0 (3) 



or (Ledoux (1947); Epstein (1979)) 



a P \ (a} + y(p.) (%) >0 , (4) 



dS) p{Yi] \dPj„ Y \ aY >Je,s\iPJ„ 

where, in the usual thermodynamic notation, the variables in the subscripts are to be held constant 
when taking the partial derivatives. If the sum involving the composition gradient can be neglected 
(as usually the case in non-degenerate stars such as the Sun and other main-sequence or giant stars) , 
the condition reduces to the familiar Schwarzschild condition (Schwarzschild (1906)) of decreasing 
entropy per baryon with increasing pressure. 



3. Thermodynamic proof 

In what follows, it is assumed that the fluid is isentropic (i.e., S = Sq = constant through- 
out) and that in the unperturbed configuration the matter is fully catalyzed, i.e., the param- 
eters {Yi(P, So)} are those corresponding to full chemical (and thermodynamic) equilibrium at 
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the local pressure. This equilibrium state is determined by minimizing the enthalpy per baryon, 
h(P, Sq, {Yi}) over all values of {Yi} allowed by conservation laws, 

h eq (P,S ) = mmh(P,S ,{Yi}). (5) 

At given entropy per baryon So, the resulting function h eq (P, So) defines a curve on the pressure- 
enthalpy plane (see Fig. 1). Point / on the curve represents the initial (equilibrium) state of the 
fluid element to be perturbed. The perturbation, at constant S = So and {Yi}, takes the fluid to 
a different pressure and out of equilibrium, therefore to a point above the equilibrium curve. The 
trajectory of the fluid element is tangent to the equilibrium curve at the point representing the 
equilibrium state of the fluid element, 



dh\ _ ( dh\ 
and curves above the equilibrium line to both sides of this point, 



(6) 



f d 2 h\ _ fd 2 h\ _ ^ / d 2 h \ fdYj\ (9Yj\ 

where the identity can be obtained from standard thermodynamic relations, together with the 
minimum condition on h in equilibrium. The inequality becomes an equality only in the special 
cases when either the equilibrium concentrations Yi do not depend on pressure, or the enthalpy is 
insensitive to changes in these concentrations. 

We note that v = (dh/dP) s ^ Yi } is the volume per baryon, so equation (7) can be interpreted 
as a condition on the derivatives of v with respect to pressure: 

dv\ ( dv \ 

In the non-relativistic limit, the mass density p = m/v, where m is the (constant) baryon mass, so 
equation (8) is equivalent to the condition for convective stability or marginal stability (eq. 3) for 
an isentropic equilibrium state. In the more general case of a relativistic fluid, the gravitating mass 
is the enthalpy (including rest mass; e.g., Misner, Thorne, & Wheeler (1973)), rather than the rest 
mass alone. Therefore p should be interpreted as the enthalpy per unit volume, p = h/v. Since, to 
lowest order, the changes in h are the same in the perturbations and in the equilibrium state (eq. 
6), the differences in p are still dominated by the differences in v, and the argument goes through 
as before. 

We conclude that isentropic, fully catalyzed matter in a gravitational field is never unsta- 
ble to convection (and usually stably stratified). Of course, the stability is enhanced if there 
is an additional, positive contribution from the entropy gradient to the Brunt- Vaisala frequency 
{{dS/dP) eq < 0), as is the case in an isothermal fluid. 
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4. Discussion 

An indirect argument leading to the same conclusion could have been given by reductio ad 
absurdum. Assume that the matter is isentropic, in chemical equilibrium, but unstable to convection. 
Then, convective overturn interchanges fluid elements initially at different pressures, which at their 
new positions are out of chemical equilibrium. Over a long enough time, reactions reestablish 
chemical equilibrium, and therefore the initial, unstable state. This process leads to a perpetuum 
mobile, inconsistent with energy conservation. This argument, although strongly arguing for the 
previous conclusion, gives no insight into the local physics that prevents the described process from 
occuring, but instead provides an additional motivation for the analysis given above. 

Let us now consider the application of the result to matter inside neutron stars. Immediately 
after the initial collapse, these stars are hot and opaque to neutrinos (Sato (1975); Arnett (1977)), 
whose non-uniform distribution violates the condition of full chemical equilibrium. Both this effect 
and a significant entropy gradient act towards destabilizing the star (Epstein (1979)), leading 
to vigorous convection (Keil, Janka, & Mueller (1996); Pons et al. (1999)). However, neutron 
stars older than a few seconds have lost essentially all their initial entropy and no longer have 
any neutrinos trapped inside (e.g., Pons et al. (1999)). At this time, they have settled into an 
essentially zero-entropy, full chemical equilibrium state. The result derived in the previous section 
implies that, regardless of the large uncertainties in the equation of state for dense matter (Lattimer 
&; Prakash (2000)), these stars are stably stratified. Marginal stability is in principle possible, but 
unlikely to hold over extended regions. Furthermore, the small remaining entropy is distributed so 
as to give a nearly uniform temperature (Nomoto & Tsuruta (1987)), which also has a stabilizing 
effect. The implications of stable stratification for different aspects of neutron star astrophysics 
have been mentioned in the Introduction. 

The writing of this Letter was motivated by a talk given by C. Thompson within the program 
on Spin and Magnetism in Young Neutron Stars at the Institute for Theoretical Physics in Santa 
Barbara. Parts of it are based on discussions with P. Goldreich which took place several years 
before. Thanks are also due to J. Veliz for preparing the figure. This work was supported by 
FONDECYT (Chile) under grant no. 8970009 (Lineas Complementarias), and by the National 
Science Foundation (USA) under grant no. PHY94-07194. 
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Fig. 1. — Schematic representation of the enthalpy per baryon as a function of fluid pressure. The 
long, curved line represents the equilibrium states (which minimize the enthalpy per baryon at given 
pressure and entropy per baryon). Points above this line correspond to non-equilibrium states, and 
there are no physical states below the line. The shorter line with arrows represents the oscillating 
trajectory of a fluid element slightly perturbed from its equilibrium position at point /. 



